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C^ We develop an analytic formalism that allows one to quantify the stability prop- 

QQ erties of X-type and Y-type junctions in domain wall networks in two dimensions. 

A similar approach might be applicable to more general defect systems involving 

^ junctions that appear in a range of physical situations, for example, in the context 

of F- and D-type strings in string theory We apply this formalism to a particular 

field theory. Carter's pentavac model, where the strength of the symmetry breaking 



CN is governed by the parameter lei < 1. We find that for low values of the symmetry 
> 

^V^ breaking parameter X-type junctions will be stable, whereas for higher values an X- 

\ type junction will separate into two Y-type junctions. The critical angle separating 



r^ the two regimes is given by a^ — 293°^y/]q and this is confirmed using simple nu- 

^"^ merical experiments. We go on to simulate the pentavac model from random initial 

^ conditions and we find that the fraction of X-type junctions to Y-type junctions is 

rS higher for smaller e, although X-type junctions do not appear to survive to late time. 

We also find that for small e the evolution of the number of domain walls A'dw in 

Minkowski space does not follow the standard oc t~^ scaling law with the deviation 

from the standard lore being more pronounced as e is decreased. The presence of 

dissipation appears to restore the t~^ lore. 
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I. INTRODUCTION 

Understanding the dynamics of junctions in networks of topological defects is of great 
interest in a large variety of physical systems. The relevance of junctions in cosmic string 
networks has recently been motivated by developments in string theory where the interaction 
of F- and D-strings [Ij can be modeled by junctions between composite (p, q) superstrings 
[SHlj . They are also relevant in condensed matter systems [5j and in building field theoretic 
models of dark energy [6]-[8]. These dark energy models are postulated to be formed by a 
frozen network of domain walls, but it has remained a challenge to build models in which 
the freezing of the network is natural. The differences between systems which favour Y-type 
and X-type junctions have been highlighted, for example [91412]. with the perception being 
that models with X-type junctions are more likely to lead to frustration. In addition it 
has been suggested to introduce a field coupled to the domain wall forming field with an 
unbroken continuous U{1) symmetry and hence a conserved Noether charge [131 E] oi* ^ 
discrete topological charge in an SU{5) x Z2 theory p^. It has been shown that both such 
models can possibly lead to a frozen network. 

In this paper we aim to develop an understanding of the stability properties of X-type 
junctions in the pentavac field theory that was first proposed by Carter [HI [17] . The 
Lagrangian is that for two scalar fields $ = l^le^"^,^ = |^|e^^ which is invariant under 
global U{1) X U{1) transformations except for the symmetry breaking term in the potential 
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COS (20 + ip) + cos {2ilj 



(1.1) 



Carter showed that this model prefers X-type over Y-type junctions in the limit e ^ 0; in 
this work we will consider a range of values of e 7^ 0. 

In recent work Avelino et al [18j studied this model and make the claim that, although 
X-type junctions do form, their eflFect is not sufficient to prevent the relaxation of the system 
in to a scaling regime where the number of domain walls, A^dw(^)5 scales like t~^. Moreover, 
they suggest that the X-type junctions which are seen are as a result of the algorithm [19] 
they employ to model the expansion of the Universe. However their simulations probe a 
limited range, taking 0.05 < e < 0.2. 

We will present a perturbative calculation of the energy associated with domains walls 
in this model which we use determine the stability of X-type junctions with particular 



intersection angle a. This allows us to determine a^^ the angle above which the X-type 
junction is stable and below which it splits into two Y-type junctions, as a function of e. 
We check this result against that found in numerical simulations of the full field theory and 
find excellent quantitative agreement. This result suggests that, for sufficiently low e, X- 
type junctions will be stable in all but the most extreme circumstances, for example, when 
junctions which have very small a are generated by the dynamics. We investigate this using 
simulations starting with random initial conditions. 

The structure of this paper is as follows. We explain what junctions are and present 
the general model we use to investigate junction stability in section |II A[ In section |IIB 



we present our perturbative analytical method we use to understand the stability of an X- 
type junction against decay into two Y-type junctions. We present Carter's pentavac model 



in section III, and simulate X-type junctions within the model. We present results from 



numerical simulations of Carter's pentavac model from random initial conditions in section 



IV and some concluding remarks in section W\ 



II. BASIC PICTURE 



A. Walls and junctions 



When a field theory has multiple disconnected vacua, the interpolation of the field between 
spatially adjacent vacuum states is a domain wall (see, for example, [20j). If more than two 
vacuum states meet at a point in space the intersection of the set of domain walls produce 
junctions. Four vacua cycled round a point produce X-type junctions, and three vauca 
produce Y-type junctions. See Figure [T] for a schematic representation of X-type and Y-type 
junctions. The tension of the walls is given by the "distance" (calculated via non-trivial 
integration) between points in the field configuration space. 

An X-type junction can decay into two Y-type junctions if the energy of the configuration 
is reduced, as depicted in Figure [l](b) . As the X-type junction splits apart a high tension 
wall Til is created, with all other walls being low tension Tj. One can think of an X-type 
junction as being two Y-type junctions "glued" together and separated by a high tension 
wall. Hence we can investigate their stability by considering whether two Y-type junctions, 
separated by a high tension wall, are in equilibrium. By simple inspection of the middle 
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(a) Y-type and X-type junctions 







(b) X-type junction decay 

Figure 1: Schematics of Y-type and X-type junctions. In (a) we show three vacua ABC cycled 
around a point to produce Y-type junction and four vacua ABCD an X-type junction. In (b) we 
give a schematic view of an X-type junction decay to two Y-type junctions. The relative thickness 
of the lines denote the relative tensions of the corresponding domain wall: the thick horizontal wall 
is a high tension wall Tn and the thin lines are low tension walls Tj. 

panel of Figure [l](b) and resolving forces, X-type junctions are stable if 



Til > 2cos(|a)Ti. 



(2.1) 



If this condition is satisfied, the horizontal high tension wall is strong enough to overcome 
repulsion of the two low tension walls, and the X-type junction retains its form. There will 
be a critical angle a^^ dependent upon the model parameter, above which X-type junctions 
will be stable and below which there will be separation into two Y-type junctions. 

B. General formalism 



We will consider a static domain wall configuration, where the field (j) (say) interpolates 
from (/)(x = +oc) = 74to(/)(x = — oc) = 5, giving rise to a domain wall between the vacua 
A and B. Integrating the energy density of the relevant field theory over the domain wall 
(between x = ±oc) gives the tension of the wall separating A and B. An alternative method 
of finding the tension is to integrate through the field manifold itself, rather than through 
the spacetime manifold. We will show how to convert between these two approaches. In 
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(a) Domain wall in spacetime 
manifold. 




(b) Domain wall trajectory in field manifold. 

Field values that minimize the potential are 

displayed by full circles. 



Figure 2: Schematic of a domain wall configuration. In (a) we show a configuration where there 
are regions of space occupying different vacua, separated by a domain wall whose energy density is 
depicted by the dotted line. In (b) we show a field trajectory between two different vacuum states, 
signifying the presence of a domain wall. 

Figure [2] we schematically depict a domain wall in the spacetime manifold v^ith coordinates 
{x^} and field manifold v^ith coordinates {(/>i}. 

Integrating the energy density <? of a one-dimensional static manifold v^ith Cartesian 
coordinate x G (— oc, +oc), having Riemannian kinetic metric G^j, gives the tension 



T = E 



dx£ 



1 dip'dip^ 



(2.2) 



where the {<y9*} are field coordinates and V is a potential function which has zero value at 
its minimum, Kain = 0. Introducing an affine parameter A in the field manifold, one finds 
that the tension can be written as 



1 



T = / dA — G,,(^'V'^ + Vx'], 



2x 



where 



dx 



X = 



, dip 

<p = 



(2.3) 



(2.4) 



dA' ^ dA' 
To eliminate the dependance of the tension on the spatial coordinate x, we minimize the 



integrand with respect to x' (i.e. setting 5T/5x' = 0), giving 



X = 



Gijif'iifi'^ 



2V 



(2.5) 



which can be substituted into (2.3) to give 



dA^/2VG~^^^^. 



(2.6) 



This expression will form the starting point of our subsequent discussions. If the endpoints 
of integration are taken to be different vacuum states, then T is the tension of the domain 
wall between the vacua. 

So far the computation of the tension has been performed over some continuous manifold, 
such as the circle S^ or torus S^xS^^ having continuous symmetry. The vacuum manifold of a 
field theory such as the pentavac model is found by constraining motion to some submanifold 
via a small symmetry breaking parameter ^, say. We now perturbatively expand the tension 



(2.6), by writing the potential and metric in orders of ^, 



V = ^y (1) + ^ V(2) , Gij = Sij + gG\f , 



(2.7) 



where 6ij is the Kronecker-delta symbol and we used the fact that 1/min = 0. Substituting 



the expansions (2.7) into (2.6), the tension can be written perturbatively as 



T ^JdAV2gV^[l + f (^ + g!]V'V'^) + Oig')^ , (2.8) 

where we have used Sij(p''^(p'^ = 1. We require V^^^g > 0, a necessary condition for this 
analytic method to work. 



C. Specific form of the energy density 

We now derive an expression that can be used to calculate the tension of a domain wall 
in models with energy density given by 
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2 iN2' 



£= [|V$|2 + |V*p]+- (|$|2_i)^+(|^|2_i)^ +^|$n*|2F((/),^), (2.9) 



where 



$ = |$|e"^, * = |*|e"^. 



(2.10) 



We now perturb about the solution for the ^ = model and integrate along the geodesic 
trajectories in the vacuum manifold to give the tension. In particular we write 



$ = (1 + g-fi) e'^, ^ = (1 + ^72) e 



ii/> 



Substituting this perturbation into the energy density (2.9) whilst keeping terms up to (9 (p) 



gives 
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(V0)2 (1 + 2^71 + q'^I) + {V^f (1 + 2^72 + o'll) + q" [(V71)' + (V72)'] 



+^F((/>, i,) + q" [2(71 + i2)F{<t), ij) + A(7i' + 72')] + 0{g'). 
This can be written as 
S 



(2.11) 



1 ^ dip' d(f^ 1 ^ dY d7-^ 

2 dx dx 2 dx dx 



^' = {0,^}, 



if one makes the identifications in terms of the expansion (2.7) 



y(i) = F(0,^), y(2) = 2(7i+72)W,^) + A(7i' + 72'), G^f = 2diag(7i,72). 



Mi. 



q\ 



i]i 



^3 



Then using ( |2.8[ ) one finds 



/"dA^/S 



2qF ll^Q 



71 (1 + 0'^) + 72 (1 + r) + ^ (7? + 72^) 



+ 0(£)2) .(2.12) 



The trajectories 7j that extremise this functional are found to be 

7i = --(l + 0'2), ^^^ 



^(1+0'^), i. = -^{i+r). 



and, using the condition 0'^ + -0'^ = 1, ( |2.12[ ) becomes 
This can be written as 



A 



(2.13) 



(2.14) 



where the zeroth and first order terms in the tension are given by 

U= fdA y^F^/\ W = A I dk y2^F^/^ ^ = - (4 + 0'^ + ^'^) . (2.15) 

To finally obtain stability criteria, one must integrate these expressions for a particular 
theory (i.e. choice of F((/), ?/;)), over particular trajectories (such as high and low tension 
walls). This is a task we leave for the next section, where we continue with a specific example 
of a field theory. 



III. SPECIFIC EXAMPLE: THE PENTAVAC MODEL 

In this section we will consider a specific field theory: Carter's pentavac model [HI [T7] . 
We will begin with a discussion of the salient aspects of the pentavac model and continue to 
give an understanding of the stability of X-type and Y-type junctions. Specifically we will 



use the perturbative expression for the tension (2.15) to obtain an analytic dependance of 



the critical internal intersection angle o^c, as a function of e. 

A. Pentavac potential 

Carter's pentavac model [161 [El has two scalar fields ^(x^), $(x^) G C whose dynamics 
are described the Lagrangian density 

1^ = \d^^d^^ + T.d^'^d^^ - V ($, ^) , (3.1) 

Zj Zj 

interacting in the potential 

y = ^ (|$|2 - r^^Y + ^ (|^|2 _ ri'Y + ^|$n^|2 (cos^ + cosx + H) , (3.2) 

where 

$=|$|e^^, ^=|^|e^^, e = 2(t) + ^, X = 2^-0, (3.3) 

and Vo is chosen so that V^in = 0, to be in accord with the requirements of the analytic 
calculation we will perform. One can easily show that by rescaling the fields according to 
$ h^ $7y, ^ h^ ^jT] and the space-time coordinates x ^ x/yXrj^ there is only one meaningful 
combination of model parameters, namely e = q/\. Hence, without loss of generality we 
will set the model parameters A = ry = 1 throughout; thus, the limit e ^ corresponds to 



q/\ -^ which can be achieved by A ^ oc. With the rescaling the potential (3.2) becomes 



V =-{\^^ -if + -{\^^ -if + e|$n^|2 (cos^ + cosx + H) . (3.4) 

So that the e < sector can be understood on the same analytical footing as the e > 
calculation presented in the previous section, we must be able to identify a transformation 
that relates the e > and e < sectors. For the pentavac model let us write 

J(6^,X, e, Vb) = 6(cos6^ + cosx + Vb), (3.5) 



then, the transformation (and invariance) is 

J{e, X, e, Vo) ^ J>, X, e, H) = J(^ + TT, X + tt, -e, -K)) = J{0, X, e, H)- (3.6) 

If e = the theory has continuous U{1) x U{1) symmetry, meaning that the vacuum 
manifold is a torus, S^ x S^. When e 7^ the continuous symmetry is broken and there are 
discrete minima on the torus. One can easily find that in the vacuum the field moduli are 

|$P = l^p = ^---. (3.7) 

' ' ' ' 1-4|6| 

This specific choice of symmetry breaking potential is periodic, and therefore it is nat- 
ural to look for transformations of the phases (0, V^) that leave the Lagrangian invariant. 
Transformations of the form 

01 — > ^ = (j) + 5^, ip \ — > ijj = 1/; + 5^, (3.8) 

leave the symmetry breaking potential unchanged if the argument of the cosine-terms in the 
potential increment by an integer multiple of 27r, i.e. if 

20 + ?/; — y2(j) + ^Jj + 2^n, 2?/; - — > 2^Jj - cj) + 27rm, (3.9) 

where n^m ^ Z and hence 

271" 27r 

5^ = — (2n - m) , 5^ = ^ (2m + n) . (3.10) 

Therefore, we can introduce two winding numbers (n, m) which can be used to construct 
transformations of the fields which leave the Lagrangian density invariant, and therefore 
construct the phase combinations that constitute the vacuum manifold. In the e < sector 
the five vacua of the theory are given in terms of the phases (0, ip) by 

(0, 0), A: (-27r/5, 47r/5), B : (47r/5, 27r/5), C : (27r/5, -47r/5), D : (-47r/5, -27r/5), 

and when e > 0, by the phase combinations 

(tt, tt), A: (-3^5, 7r/5) B : (^5, 3^5), C : (3^5, - ^5), D : (-7r/5, -37r/5). 



We present an illustration of the vacuum manifold in Figure p^ in the e > sector (the 
e < manifold has the same structure, but with the maxima and minima reversed). The 
plot we present is in the space of field phases, (0, '0), where the figure is periodic in both 
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Figure 3: The vacuum manifold of the pentavac model in the e > sector. We have marked on two 
wall trajectories between the vacua B and C; the low tension interpolation, Tj, and high tension, 
Tii, (remembering that this grid has the topology of a torus). The minima of the theory correspond 
to the blue (dark) regions, and the maxima to the orange (light) regions. 

directions (i.e. toroidal). Marked onto the figure are two energetically diflFerent trajectories 
between the vacua B and (7, i.e. the two different types of domain wall: high tension Tu 
and low tension Tj. Two such trajectories exist between all pairs of vacua. 

In the limit e ^ the vacuum |$| = |^| ^ 1, which has the interesting consequence of 
decoupling the field phases; that is, taking |$| = |^| = 1 the Lagrangian density can be 
written to leading order in e as 

£ = ^(1 - 4e) [d.Od^O + d.xd'x] " 2e {cos^{9/2) + cos2(x/2) - l) + 0{e^). (3.11) 

This describes two non-interacting sine-Gordon kinks and provides us with a very useful 
result: in the asymptotic e ^ limit an X-type junction is absolutely stable. 

B. Analytic calculation of ac{e) 



If we define / = J/e then for e > 0, /((/), t/^) = cos0 + cosx + 2 (n.b. for e < 0, 



f{(j)^i/j) = COS0 + cosx — 2) and we used this to compute (2.15). Now we pick typical high 
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and low tension trajectories and perform the necessary steps to compute ac{e): 

• Low tension Tj: the low tension trajectory has x = ^r, so that fi = 2cos^(0/2) and 
dA = -y^dO. The variation of the field phases 0, ip along a low tension trajectory is 
found by inspecting the vacuum manifold and deducing the equation of a line linking 
vacua. One finds 0'^ = 1/5, V^'^ = 4/5. 

• High tension Tu: the high tension has x = ^^ so that fu = 4cos^(6^/2) and dA = 
-^ V^d0. The variation of the field phases along a high tension trajectory is found to 
be 0^2 _ 1/10, 7/;'2_ 9/10. 

From these expressions one can find that the zeroth order tensions satisfy 



Uu = 2t/i. 



(3.12) 



Using these ingredients in (2.15), one can easily compute 



Wn 1928 



Wi _ 468 
T^ ~ 150' 



(3.13) 



f/ii 300 ' 
The stability of an X-type junction is assured if the internal intersection angle a satisfies 



(2.1). We take ttc to be the intersection angle at equality. In the notation of (2.15), we can 



expand to first order to find 



T, 



II n \ Uu 

—— = cos iftc = —rr 
2Ti ^2 c; 2C/j 



Wn _ Wi 



Using our computed values (3.13), this becomes 



248 
cos (l^c) = 1 - e— - + 0{e^). 



and expanding cos {^a^^ '^ 1 — \a^^ one finds 



0{e') 



(3.14) 



«c ~ 293°^, 



(3.15) 



which holds for < e <C 1. In fact, we see that by the transformation/invariance relationship 



(3.6) we can extend (3.15) to include e < 0, and write a^ ^ 293°A/|e| . Hence, we see that 



the expression predicts that systems with low e are more stable to acute perturbations that 
those systems with high e. 
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C. Numerical simulation of X-type junctions 

In order to check the veracity of the analytic calculation, we have performed a suite 
of simulations evolving the equations of motion. We will investigate the stability of X- 
type junctions by constructing an X-type junction with a particular a and investigate the 
conditions under which it decays, if indeed it does decay 

In our numerical investigations we discretize space to fourth order, and time to second 
order - evolution is performed with a leapfrog algorithm, on a grid of A^^; x A^^ grid-points. 
We use space step-size Ax = 0.5 and time step-size At = 0.1. To smooth out initial 
unphysical discontinuities we damp the fields for the first 200 time-steps (i.e. t < 20) and 
then remove damping, allowing the system to freely evolve. All simulations use periodic 
boundary conditions. We will display images of the field, where each colour corresponds to 
the vacua closest to the field at a given location. The initial conditions for our numerical 
simulations have X-type junctions with a "tunable" initial internal angle a by choosing 
a = 2tan~-^ (Ny/Nx) as shown in Figure El 

To make the model numerically tractable we cast each complex field as two real fields, 
$ = 01 + 102, ^ = 03 + 104, and then one finds that the potential can be written as 

-'- ^ 2 , /2 i\2 , 1 / /2 , /2 i\2 



V{<|>^) = ^(0? + 0i-l) +-{<PI + Cf>l-1) 



+€ 



(3.16) 



In Figure [5] we present the evolution of the field for an initial configuration with an X- 
type junction having an initial internal angle a = 53° (achieved by setting the number of 
horizontal and vertical grid-points to A^^; = 256, Ny = 128, respectively). The case of e = 0.1 
is in the top row. The X-type junction clearly breaks apart into two Y-type junctions. 
However, the e = 0.03 case presented in the bottom row, the X-type junction retains its 
form. Clearly ac = 53° for some value of e between e = 0.1 and e = 0.03. 

A plot of the distribution of the phases (0, ip) at a given time will show where the points 
are located on the vacuum manifold, and will highlight any fundamental differences between 
systems whose X-type junctions have and have not decayed. Referring back to the vacuum 
manifold in Figure [3j we note that if any phases are found across the origin then high tension 
walls exist. Conversely, if this is not the case no high tension walls exist between opposing 
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Figure 4: The initial field configuration for setting up the X-type junction under investigation. 
There are 8 equivalent X-type junctions in this configuration - with the figure being periodic in 
both directions. 




^S^^^S 




Figure 5: Evolution of the field in the pentavac model, at times t = 0, 10, 20, 30 (left to right) for 
e = 0.1, 0.03 (top to bottom) and initial internal angle a — 53°. Each colour represents one of the 
vacua in the theory. For e = 0.1 the X-type junction with a — 53° is not the lowest energy state 
and it breaks apart into two Y-type junctions. In the e = 0.03 case the X-type junction is the 
energetically preferred state and the X-type junction retains its form. 

vacua and adjacent vacua must wrap the torus to produce high tension walls. We present 
the phases for e = 0.1 and e = 0.01 in Figure |6} for an initial X-type junction with a = 53°. 
One can clearly observe that there exists high tension walls in e = 0.01 case, and only low 
tension walls in the e = 0.1 case. These images represent the final state of the system. 

The potential energy density isosurfaces for the case of e = 0.1 are presented in Figure 
[7| The value of the potential energy density is denoted by the relative brightness at a given 
location. We see that in the second image a high tension wall has formed, which subsequently 
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V^ 




ip ^' 




Figure 6: The phases of the fields, (0, '0), evaluated at each point in space when t = 30 following 
relaxation of a configuration having an initial X-type junction with a — 53° for two different values 
of the symmetry breaking parameter e. The case of e = 0.01 has retained its form as an X-type 
junction, while the case of e = 0.1 has relaxed into a Y-type junction configuration. One can clearly 
see that one crosses the origin in the e = 0.01 case, which correspond to high tension walls, whereas 
this does not take place in the 6 = 0.1 case with all walls being low tension. 




Figure 7: Images of the potential energy density at times t = 0, 5, 10, 15 (left to right) for a decaying 
e = 0.1 X-type junction. The value of the potential energy at a given location is denoted by the 
relative brightness of the colouring. One observes a high tension wall initially being generated, but 
this swiftly decays into a low tension wall. 



decays to a low tension wall. This is compatible with the schematic decay process depicted 
in Figure [l](b) . At first sight it appears that there is some sort of discontinuous jump as an 
X-type decays into a Y-type junction: the field trajectories transition from wrapping around 
the toroidal vacuum to not wrapping. This transition is energetically driven, and is therefore 
not a cause for concern. 

We have developed an automated algorithm to compute the critical angle, a^^ at which 
an X-type junction becomes unstable as a function of e. We first compute the length and 
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Figure 8: The values of ac (in degrees) deduced numerically for various symmetry breaking pa- 
rameters e in the pentavac model. The solid hne is the best fit to the data, whereas the dotted line 
is that predicted by the theory presented in the first section. 



width of the wall for the junction at the centre of the grid, L and S respectively. If the 
junction is of X-type then the walls will be of comparable length, L/6 ^ 1; however, if the 
junction is of Y-type then the horizontal wall will be much longer than the vertical wall, 
L/6 ^ 1. We compute the value of L/6 after simulation time t = r^ for initial conditions 
with an internal intersection angle a. If L/6 < C, the junction is deemed to be of X-type 
and the vertical size of the box is reduced to lower a. This process continued until a Y-type 
junction with L/6 > ( is found. We take r = 30 for A^^^ = 256, and ( = 2.5 whose value was 
determined by convergence tests. 

We present numerical estimates of ac in Figure [8] for a range of positive and negative 
values of the symmetry breaking parameter. Our algorithm suffers from resolution effects as 
e ^ since the domain wall width widens, and therefore we only show results for \e\ > 0.02. 
The results show that: (1) ac is largely independent of the sign of e, something which is 
understood by the structure of the vacuum manifold (the manifolds for e > and e < 
have identical structure); (2) the numerical results have the same dependance on e as the 
analytical prediction, however, their values are ^ 10% lower which can be understood as a 
rather trivial numerical resolution effect - the calculation predicts the value of ac for the 
moment an X-type junction breaks apart by an infinitesimal distance whereas the numerical 
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calculation has to wait until L/5 = 2.5. Unfortunately we are unable to reduce C, any more 
due to numerical resolution effects. 

We note that it is not possible within our framework to perform the reverse process. It 
is not energetically possible to start off with two Y-type junctions and for them to contract 
into a high tension wall. 

IV. SCALING DYNAMICS FROM RANDOM INITIAL CONDITIONS 

In the previous section we have shown that both X-type and Y-type junctions can be 
stable, dependant on the intersection angle a and e. We now develop this further by testing 
how the dynamics of random domain wall networks are affected by varying the symmetry 
breaking parameter e. From our investigation of idealized X-type junctions in the previous 
section we have seen that a lower e implies that an X-type junction is stable for a wider 
range of intersection angles, for example if e = 0.01 then X-type junctions will be stable for 
a > 30°. Performing simulations of a random network will allow us to put a measure on the 
probability distribution of a in realistic situations. 

We have evolved the equations of motion in Minkowski space with random initial domain 
occupation on lattices with P = 4096 grid-points in each direction, and space step-size 
Ax = 0.5 and an initial correlation length of one grid-square. We apply a damping term 
of magnitude 0.5 for the first 200 time-steps to smooth out unphysical discontinuities that 
result from our initial conditions. We present images of the resulting field configurations as a 



function of time in Figure |9] and Figure [T0| for various values of e. The "mottled" appearance 
of the e < 0.01 simulations is due to the rather shallow potential wells defining the domains 
as e ^ 0; however domains still clearly form. By carefully inspecting the field images we 
have been able to find that there is a higher fraction of X-type junctions to Y-type junctions 
in simulations with lower e; however, the X-type junctions that are present at early time do 
not seem to survive to late time. 

To understand how the wall network evolves we calculate the number of domain walls at 
every time-step and build up a plot of the evolution of the number of domain walls N^^{t) 
and compute the scaling exponents 7 defined by A^dw oc t^. A domain wall is numerically 
defined when the vacuum that the field occupies at a given lattice site is different from the 
vacuum occupied by any of its neighboring lattice sites. The standard lore is 7 = — 1 [^T-'^S] 
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Figure 9: Evolution of fields in the pentavac model, for e = 0.1, 0.05, 0.03, 0.01 (top to bottom) from 
random initial conditions on a P = 4096 square grid. Images are at t = 80, 160, 320, 640, 1280 (left 
to right) and the colours correspond to the vacua of the pentavac theory. Each colour represents 
one of the vacua in the theory, and are assigned by finding which of the vacua the field is closest 
to at each location. 



in which case the network collapses as fast as causality allows. Deviation from this value 
signals the existence of some force resisting the collapse of the wall network and the possible 
formation of a lattice. 

We plot the evolution of the number of domain walls for the systems in Figure |9] and 



Figure 10 in Figure fllf a). For e > 0.03 the number of domain walls follows the standard 



t ^ power law. However there is substantial deviation for e < 0.01. Plotting the evolution 
of the scaling exponents j{t) in Figure [TT](b) confirms this. 
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Figure 10: Evolution of the fields in the pentavac model, for systems with e = 
0.01,0.008,0.005,0.001 (top to bottom) at times t = 80,160,320,640,1280 (left to right), on a 
grid with P = 4096 lattice sites and space step-size Ax = 1.0. 



One might be concerned that this effect is due to the finite resolution and dynamic 
range of the simulations particularly since the width of the domain wall is expected to be 
S^ oc 6~2 . To test our results we have performed a number of convergence tests. Simulations 
with various space step sizes, correlation lengths and a variety of different physical sizes of 
grids were used. A correlation length ^ is imposed by making the system have blocks of ^x^ 
grid-points initially in the same vacuum. We will focus on the case of e = 0.01 since it is close 
to the parameter range where the non-standard behavior starts to become noticeable. First 
we evolved with P = 512, 1024, 2048 with correlation lengths ^ = 1, 2, 4 and space step-sizes 
Ax = 0.5,0.25,0.125 respectively; these fix the physical size of the box while varying its 
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(a) Evolution of the number of domain walls, A^dw(^)- 




e = 0.01 - 
8 = 0.03 
8 = 0.05 

8 = 0.1 



600 
Time 




(b) Evolution of the scaling exponents, j{t). 

Figure 11: Evolution of the number of domain walls A/^w and scaling exponents 7 from random 
initial conditions, in the pentavac model, for various symmetry breaking parameters e, where the 
grid size is P = 4096. In the left-hand panels. Ax = 0.5 is used, with scaling exponents being 
computed in bins of 100 units of time. In the right-hand panels. Ax = 1.0 is used and scaling 
exponents are computed in time-bins of 200. One can easily observe that the scaling dynamic 
N^w{t) is different in the cases e > 0.01 and e < 0.01. The exponent, 7(t), of the power-law that 
^dw(^) satisfies clearly changes as e changes, with 7 ^ as e ^ 0. 

resolution. Secondly we change the physical size but retain the number of grid-points. To 
ensure that the Courant condition is satisfied we changed the time step-size for systems with 



Ax = 0.125 to At = 0.05. In Figure 12 we show a plot of the evolution of the number of 
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(a) Fixed physical size 



(b) Fixed P = 2048 



Figure 12: Evolution of the number of domain walls with e = 0.01, with varying resolution: (a) 
P = 512, 1024, 2048 with varying space step-size Ax = 0.5, 0.25, 0.125 and initial correlation lengths 
^ = 1, 2, 4 respectively; (b) P = 2048 and varied Ax, keeping ^ = 1. 

domain walls for these test systems. All are consistent with the same deviation from 7 = — 1 
shown in Figure [9[ and all have qualitatively the same scaling dynamics. 

A separate concern is that the effect is due to the energy that was originally localized on 
the wall network as walls collapse. As previously pointed out the energy required to created 
a wall is very low for small e. As a wall network collapses energy becomes released from the 
walls into radiation which may be sufficient to create new "secondary" walls. As the original 
wall network collapses the amount of radiation increases with a consequent increase in the 
rate of "secondary" wall production. Such a mechanism could account for the deviation 
from the standard scaling behavior we observed. 

If this "background radiation" can be removed by some means then the scaling dynamics 
may revert to the standard t~^ law. This has been tested by evolving a modified set of field 
equations, where radiation is removed. This can be done by including a constant damping 
term (but very small) or by using a modified version of the algorithm by Press, Ryden and 
Spergel (PRS), where the damping term is time-dependant, which has been suggested as a 
way of modeling an expanding Universe [HI [I9]. These equations of motion are of the form 

(4.1) 



V<j)-V^(/)+^=0. 
dcp 



In each of the modified evolution algorithms the coefficient V of the damping term takes on 
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(a) e = 0.1 



(b) e = 0.01 



Figure 13: Evolution of the number of domain walls for three "different" evolution algorithms with 
e = 0.1, 0.01; we use Ax = 0.5, At = 0.1, P = 4096 and random initial conditions. The line denoted 
"E-L M" is the evolution of the standard Euler-Lagrange equations in Minkowski spacetime. The 
line denoted "constant" has a constant damping. The hne denoted "PRS" uses the PRS algorithm 



in the radiation era. The implementation of these algorithms is given in (4.2). 



the following form: 



V 



PRS 



n 



constant 



Do t < 20, 
10-3^0 t > 20. 



(4.2) 



We take Vq = 0.5. The Euler-Lagrange equations in Minkowski spacetime have V = 0.5 for 



i < 20 and V = thereafter. In Figure 13 we compare the evolution of the number of domain 
walls in these "constant damping" and PRS algorithms (as well as with the standard Euler- 
Lagrange equations in Minkowski spacetime). It is clear that in the e = 0.01 simulation 
the evolution is substantially different in these damped algorithms to the evolution in the 
Euler-Lagrange equations in Minkowski spacetime. 

Because the inclusion of a dissipative term restores the t~^ scaling law (at least for a 
substantial period of time) it is no longer obvious which evolution algorithm corresponds to 
results relevant in cosmology. This issue may be resolved by an extensive systematic study 
of the results of the equations of motion for a range of models with as large a dynamic range 
as feasible. 
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V. CONCLUSIONS 

In this paper we have furthered the understanding of the stabihty of X-type and Y- 
type junctions. We have calculated that if an X-type junction has an internal intersection 
angle a, then ii a > ac the X-type junction is stable, whereas if a < ac the X-type junction 
decays into two Y-type junctions. In Carter's pentavac model the dependance on the critical 
intersection angle o^c upon the models symmetry breaking parameter e has been derived to 
be 

ac = 293°^^. (5.1) 

The validity of this calculation has been verified by numerical experimentation, and is found 
to hold regardless of the sign of e. 

A consequence of this relationship is that if e is relatively large (for example, e = 0.05) 
the X-type junctions cannot withstand much "squashing" before they break apart; they can 
only withstand squashing to an angle oi a ^ 65°. Therefore, X-type junctions cannot be 
expected to survive in physically relevant systems. This is what was observed by Avelino 
et al [IE] but they expanded the statement to include all systems with X-type junctions. 
However, if e is taken to be small (for example e = 0.01) X-type junctions are stable for 
intersection angles above a ^ 30° and therefore could be expected to survive in physically 
relevant systems. The ability of an X-type junction to retain its form under much more 
extreme intersections is increased as e is decreased. 

Whilst X-type junctions become more frequent in simulations from an initially random 
configuration of vacua when e is small, the X-type junctions do not survive to late time. 
For systems with e < 0.01 the scaling dynamics of the resulting network becomes modified 
compared to that of a system with e > 0.01. The amount by which the evolution is modified 
was clearly shown to depend on the value of e and we believe that this effect is due to 
"background radiation" being sufficient to create walls. 
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